On the Cooling of Electrons in a Silicon Inversion Layer 
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Cooling of two-dimensional electrons in silicon-metal-oxide semiconductor field effect transistors 
is studied experimentally. Cooling to the lattice is found to be more effective than expected from the 
bulk electron-phonon coupling in silicon. The extracted heat transfer rate to phonons at low tem- 
peratures depends cubically on electron temperature, suggesting that another coupling mechanism 
(such as piezoelectric coupling, absent in bulk silicon) dominates over the deformation potential. 
According to our findings, at f 00 mK, the electrons farther than ~ 100 fim from the contacts are 
mostly cooled by phonons. Using long devices and low excitation voltage we measure electron resis- 
tivity down to electron temperature ~ 100 mK and find that some of the "metallic" curves, reported 
earlier, turn insulating below ~ 300 mK. This finding renders the definition of the proposed 2D 
metal-insulator transition questionable. Previous low temperature measurements in silicon devices 
are analyzed and thumb rules for evaluating their electron temperatures are provided. 



Since the scaling theory for non-interacting electrons 
was constructed in the late seventies [Q, it was conjec- 
tured that a realistic two-dimensional electron system is 
insulating in the sense that its resistance diverges at low 
temperatures due to quantum interference and interac- 
tion effects. As the temperature is reduced, a logarith- 
mic resistance increase due to weak localization is ex- 
pected, followed by an exponential resistance divergence 
once strong localization commences [jjj . It was later sug- 
gested Q that strong Coulomb interaction may counter- 
act localization and lead to the existence of a two dimen- 
sional metallic state and a metal-insulator transition at 
zero temperature. 

Early experiments on various types of two dimensional 
electron gases (2DEG) supported localization. For a 
comprehensive review of early theoretical and experimen- 
tal results see Refs. [||4|. It was therefore unexpected 
when the existence of a two dimensional metallic phase 
in high mobility silicon Metal Oxide Semiconductor Field 
Effect Transistors (MOSFET) was reported in 1994 [§,§. 
The ultimate signature of such a metallic phase is resis- 
tance saturation to some residual value as the tempera- 
ture approaches zero |7|,|| . Since T — is experimentally 
inaccessible, all claims for a metallic phase rely on ex- 
trapolation from finite temperatures. Two traps lurk an 
experimentalist in this procedure: (a) The assumption 
that the observed resistance saturation indeed persists 
to zero temperature may turn wrong, (b) The electron 
temperature may (and does indeed) depart from the lat- 
tice or bath temperature, Tb, as Tf, — > 0. The latter point 
is particularly acute in silicon MOSFETs due to the high 
intrinsic contact resistance and the weak electron-phonon 
coupling. 

The experiments described here focus on the second 
point. Using the sample itself as a thermometer we 
analyze electron cooling as a function of mixing cham- 



ber temperature and excitation voltage. The main find- 
ings are: (a) Electron cooling to the lattice is more 
effective than expected from the known bulk electron- 
phonon coupling in silicon. Power dissipation is given 
by aT 3 + bT 5 with a = 2.2 x 1(T 8 W/K 3 cm 2 and 
b = 5.1 x 10 -8 W/K 5 cm 2 and, hence, at temperatures 
below ~ 0.6 K, dominated by the T 3 term g]. (b) Hot 
electron diffusion (ED) to the contacts is unimportant in 
large devices. In short Si MOSFETs, according to our es- 
timates, ED provides cooling or heating depending on the 
contact resistance, (c) Comparing cooling by phonons 
with cooling by ED at e.g., 100 mK, we find the latter 
mechanism is typically important only for electrons closer 
than ~ 100 fim to the contacts. In longer devices, cool- 
ing proceeds through phonons. (d) Cooling the electrons 
down to 100 mK solely by phonons requires an excitation 
power estimated to be less than 2 x 10 _11 W/cm 2 . For 
a 1 mm square of 2DEG with 30 Kfl resistivity this fig- 
ure implies an excitation voltage below 80 /iV. A 50 /im 
square with the same resistance requires an excitation 
voltage below 4 /zV. 

Four terminal, lock- in ac measurements at 11 Hz were 
carried out on samples similar to those used in Refs. 
The samples were 5 mm long, 0.8 mm wide Hall bars with 
2.5 mm separation between potential probes. The ox- 
ide thickness was 200 nm. Battery powered electrometer 
preamplifiers were employed to minimize electromagnetic 
interference and spurious offset bias across the potential 
probes. With 300pF cable capacitance, samples up to a 
few Mf2 total resistance were measured with negligible 
"out of phase" signal. The excitation voltage across the 
whole sample was maintained below 100 /zV to guarantee 
tolerable heating (this figure is justified below). In order 
to reduce the contact resistance we applied a larger pos- 
itive front gate voltage and — 15 U substrate bias, thus 
maintaining the same carrier concentration range. The 
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substrate bias depletes both the channel and the con- 
tacts, but it's effect on the latter is smaller due to heavy 
doping in the contact regions. This way we were able to 
reduce the contact resistance below 20 Kfi for all "metal- 
lic" densities. The substrate bias results in a moderate 
reduction of the peak mobility from e.g., 3 x 10 4 cm 2 /Vs 
at zero substrate bias to 2.1 x 10 4 cm 2 /Vs (both mea- 
sured at T = 0.3 K ). The mobility reduction is attributed 
to enhanced surface scattering due to shift of the 2DEG 
towards the Si — SiC>2 interface. The substrate bias was 
changed in the presence of red LED illumination, which 
heats the crystal, thus facilitating acceptor recharging in 
the depletion layer. The high contact resistance prob- 
lem is further complicated below 1.2K by superconduc- 
tivity of the Al contact metallization. All our "zero mag- 
netic field" measurements were, therefore, carried out in 
the presence of a 0.01 T parallel magnetic field which 
quenched the Al superconductivity. 




FIG. 1. Temperature dependence of the resistivity. 
Densities from top to bottom are n — 1.20 to 1.44 (in 
0.03 steps), 1.56, 1.8 and 2.04xlO n cm 2 /Vs. Inset: Mag- 
netoresistance at n = 1.44 x 10 11 cm 2 /Vs. The weak lo- 
calization negative magnetoresistance peak and the onset 
of Shubnikov - de Haas oscillations are clearly seen. The 
peak width decreases and the oscillations amplitude in- 
creases without saturation as Tt decreases to 0.1K. 

The resistivity vs bath temperature for a 2.1 x 
10 4 cm 2 /Vs peak mobility sample at various densities 
is depicted in Fig. |l} For temperatures above 0.4 K the 
traces are similar to those reported in Refs. |^||. At 
low densities the sample is insulating, while at higher 
densities the resistance decreases as the temperature is 
reduced, a behavior previously assigned to a two dimen- 
sional metal. For n = 1.29 x 10 11 cm -2 the resistance is 
practically independent of T in a wide temperature range, 
0.4 -T- 4 K. Such a phenomenon was previously identified 
as a metal to insulator transition (MIT) Upon 
cooling the sample below 0.4K, we find a novel result; 
the resistance of some of the "metallic" curves turns up- 



ward, showing an insulating behavior in all three devices 
studied by us. 

Measurements on a lower mobility sample with 2.5 x 
10 3 cm 2 /Vs peak mobility (Fig. |) at 0.3 K yield results 
similar to those reported earlier [3 10 for the same type 
of samples. As the temperature decreases, the low mobil- 
ity sample displays a pronounced insulating behavior at 
all measured densities, above 6 x 10 n cm~ 2 . A relatively 
small resistivity drop with cooling is again followed by a 
resistivity increase at lower temperatures. Such curves 
have been observed earlier for 2D holes in Si-Ge JlT| , 
2DEG on vicinal Si-surfaces @, high mobility Si MOS- 
FETs at higher densities, n > 2 x 10 12 cuT 2 jjjfij, and 
in 2D systems on other materials We attribute 

this resistance upturn to interference and Coulomb inter- 
action effects that grow as the temperature is reduced. 
The interference is clearly manifested in the weak local- 
ization magnetoresistance curves depicted in the inset to 
Fig. 1. 

To gain insight into the electron cooling mechanisms 
we studied the resistivity dependence upon excitation. 
We applied a small ac signal superimposed on a dc bias 
and measured the differential resistivity as a function of 
this bias. The ac excitation was typically about 100 pA, 
generating less than 65 /iV ac voltage drop across the 
5 mm sample. The analysis assumes the resistivity de- 
pends on the electronic temperature alone and uses the 
sample resistivity as a thermometer. We chose to use 
the low mobility sample due to its moderate resistance 
and monotonic temperature dependence, which facilitate 
accurate ac measurements. 
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FIG. 2. Resistivity dependence upon temperature for a 
2.5 x 10 3 cm 2 /Vs peak mobility sample. Densities given 
in 10 11 cm 2 /Vs units. 

Inset (a) to Fig. ^ depicts the channel resistivity ver- 
sus dc current at various bath temperatures. The overall 
dissipated power per unit area of the 2DEG consists of 
Joule heating, P c i = J 2 x p, by the driven current and 
heating by spurious excitations, Pq, due to e.g., absorp- 



2 



tion of electromagnetic radiation. The heat is dumped to 
phonons and to the contacts via electron diffusion. We 
ignore for the moment electron diffusion - an assumption 
to be later justified for our samples. In the general case, 
electron-phonon interaction may be due to both piezo- 
electric and deformation potential coupling with their 
characteristic T 3 and T 5 dependencies , respectively. 
Similar functional form was recently found for holes in 
SiGe heterostructures Q . We therefore characterize the 
electron phonon interaction by the coupling parameters 
a and b. 



Pd 



P d + P = a(T 3 



T 3 ) + b(T* - T b 5 ) 



(1) 



We solve Eq. 1 to find the electron temperature, T e i, 
and fit a and b in such a way that all R(T e i) data in 
inset (a) to Fig. 3 collapse to a single curve. Figure || 
shows the best fit obtained with a = 2.2 x 10~ 8 W/K 3 cm 2 
and b = 5.1 x 10 _8 W/K 5 cm 2 . Pq cannot be determined 
in this scaling procedure. The piezoelectric contribu- 
tion is, hence, found to dominate electron - phonon cou- 
pling at all temperatures below 0.6 K. This effect is un- 
expected since a Si lattice possesses inversion symmetry 
and, hence, should display deformation potential cou- 
pling only. We believe the piezoelectric contribution, 
found in our experiment, originates from the Si — SiC>2 
interface which lacks inversion symmetry. 
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FIG. 3. Channel resistance scaling for the low mobil- 
ity sample. The abscissa depicts P — Pd + aT :i + bT 5 . 
The fitting parameters are: a = 2.2 x 1CT 8 W/K 3 cm 2 
and b = 5.1 X 10" 8 W/K 5 cm 2 . Colors correspond to the 
same curves as in inset a. Insets: a) resistivity versus dc 
current at different base temperatures: 90, 120, 150, 180, 
210, 240, 330, 420, 510, 600 mK. b) the electronic tem- 
perature calculated using eq. |l| as a function of dissipated 
power for T b = 0, P = 0. 

Heat transfer due to ED can be incorporated into the 
model using the Wiedemann- Franz law Jl9| . The power 
balance reads 
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1 f irki 



d 2 T 2 

dx 2 



(2) 



where a is the sample conductivity. At temperatures 
T e i <C 0.6K the deformation potential contribution can 



be ignored. Defining T 3 = T b 3 + (P d + P )/a and 
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~3fT 



\J a/aTo, 



Eq. |2j takes the dimensionless form 

3d 2 e 2 
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2 d X 2 



(3) 



(4) 



where O = T c i/Tq and \ = x/X. A sets the character- 
istic length scale on which the T e \ approaches To- The 
only temperature scale in the problem is To, the electron 
temperature far away from the contacts, where cooling 
is dominated by phonons. Equation |I] can be integrated 
exactly. Away from the contacts, where « 1, « 
l + exp(- X ). For6»l, 9( X ) « 9 C /(1 + V(p0) X ) 2 , 
where C is the normalized electron temperature near 
the contact. 

For T = 100 mK and a = 3.3 x lO^Ohm" 1 (typical 
critical region conductance), A ps 110/im. For our 5mm 
long samples, ED is, hence, negligible and Eq. [j] can be 
safely used to estimate To. The resulting electron tem- 
perature as a function of dissipated power, is plotted in 
inset (b) to Fig. |^. Assuming that the electron-phonon 
coupling does not depend significantly on impurity scat- 
tering we use the same parameters to analyze the high 
mobility sample data. Equation ^ shows that the over- 
heating due to the measurement current is negligible and 
does not exceed 1 mK for a 50 mK base temperature. 
We have used the magnetoresistance curves shown in the 
inset to Fig. 1 to estimate the lowest sample temperature 
by fitting the curves to weak localization theory. The ex- 
tracted dephasing rate scales linearly with temperature, 
in agreement with theory Q, down to the lowest ana- 
lyzed temperature, = 100 mK. Although we use weak 
localization theory for resistivities beyond its applicabil- 
ity, the linear dependence of the dephasing rate upon 
bath temperature indicates that electrons cool down to 
T el - 100 mK. 

We turn now to investigate the implications of our re- 
sults to previous measurements. Surprisingly, there are 
only a few published experiments below 300 mK. Those 
experiments were carried out on two types of samples. 
The first type, in which the MIT was originally observed 
H||, is a conventional MOSFET that, in the lack of back 
gating, suffers from high contact resistance, especially at 
low electron densities and temperatures. These samples 
are several mm long, and heat transfer by electron dif- 
fusion to the contacts can be safely neglected. In the 
interesting regime, where the resistance drop is signifi- 
cant, the contact resistance is in the Mfi range. Such a 
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high contact resistance dictates either dc or very low fre- 
quency measurements which require, due to 1// noise in 
the amplifiers and the sample, substantial excitation cur- 
rents. Overheating is evident in both RvsT and lni? vs 
1/T dependencies (see e.g., Figs. 1 and 5 in Q and Fig. 3 
of ]2(|), even in the insulating regime where it should be 
less significant. In previous experiments pl| , p2[ , where 
T e i was measured in the "metallic" regime, it was never 
lower than w 200 -r 300 mK. 




.1 1.2 1.3 1.4 

FIG. 4. The difference between the two terminal re- 
sistance, i?2, and the ideal resistance, R2 = p x L/W, 
corresponding to a uniform sample temperature and zero 
contact resistance. = 80 mK and the excitation volt- 
age < 65 pV . Note the difference changes sign at approx- 
imately the MIT density. 



In another type of samples |£3|j24|, the gates, adja- 
cent to the contacts, are separated from the gate over 
the central part of the sample by w 50 nm gaps. This 
arrangement facilitates a higher carrier concentration in 
the contact regions and, hence, a reduced contact resis- 
tance. In contrast to the millimeters long samples, em- 
ployed here and in Refs. [Em, the effective sample length 
between the gaps in Ref. 24] is only 120/im. Substituting 
A = 60 pm into Eq. || and solving for To we find that cool- 
ing is mostly provided by phonons for all temperatures 
above ~ 340 mK. Since cooling by phonons is suppressed 
only slightly faster than by ED (T 3 compared with T 2 ), 
phonons, neglected in the analysis of Ref. [E4| done in 
Ref. [pi, are important also at lower temperatures. Elec- 
tron diffusion from the contacts may either cool or warm 
electrons in the central part of the sample, depending 
on the voltage drop, U c , across the « 50 nm gap. For a 
given electron temperature near the contact, T c , one can 
estimate the voltage drop, U c , across the gap resistance, 
p c , for which the net heat flow between the contact and 
the sample vanishes. Using Wiedemann-Franz law, the 
heat flow balance reads 



(5) 



For T b = one obtains U c /T c = (7rfc h /V3e) ~ 
156/xV/K, independent of p c . For lower (higher) volt- 
ages across p Cl electrons are cooled (warmed) by the 
contacts. For 100 mK, the bias at which the contacts 
neither heat nor cool the sample equals 15.6 pV. At 
that bias sample cooling is provided by phonons. Us- 
ing Eq. 1 we find that 100 mK electron temperature in 
a 30kf2 per square sample requires less than ~ 4 pV 
voltage drop per 50 pm square, which for the sample in 
Ref. |24j ], translates into 10 pV between the gate gaps. 
A 100 mK electron temperature throughout the sample 
thus requires less than ~ 40 pV (50 mK requires < 20pV 
) bias across the whole sample (neglecting the voltage 
drop across the high density contacts). Electron over- 
heating in the contact regions occurs in our samples, as 
illustrated in Fig. ||. The contact regions are always more 
resistive than the rest of the sample in the "metallic" 
regime, n > 1.29 x 10 n cm~ 2 , and less resistive in the 
insulating one, n < 1.29 x 10 n cm~ 2 , indicating they are 
warmer than the rest of the sample. 

In summary, we found that the cooling of electrons 
by phonons in Si MOSFETs is substantially more effec- 
tive than that expected from the bulk electron-phonon 
coupling. Long samples arc cooled by phonons. Short 
samples are susceptible to heat transfer by ED from the 
contacts. Low temperature measurements in short sam- 
ples require low excitations. Using 5 mm long devices, 
where electron heating in the contacts can be neglected, 
we were able to measure electron resistivity down to 
~ 100 mK temperatures. We found that p(T)-curves in 
the well pronounced "metallic" regime, corresponding to 
(1.3 — 1.44) x 10 cm~ 2 densities and a three fold re- 
sistance reduction with cooling, turn insulating below a 
certain temperature. That temperature is ~ 300 mK at 
p ~ h/e 2 and decreases as the density is increased and 
p decreases. Our finding renders the definition the 2D 
metallic phase and the metal-insulator phase transition 
questionable. Further work is needed to determine the 
physics behind the observed resistance upturn. 
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